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Abstract 


The existence of black holes and of spacetime singularities is a fun- 
damental issue in science. Despite this, observations supporting their 
existence are scarce, and their interpretation unclear. We overview how 
strong a case for black holes has been made in the last few decades, and 
how well observations adjust to this paradigm. Unsurprisingly, we con- 
clude that observational proof for black holes is impossible to come by. 
However, just like Popper’s black swan, alternatives can be ruled out or 
confirmed to exist with a single observation. These observations are within 
reach. In the next few years and decades, we will enter the era of precision 
gravitational-wave physics with more sensitive detectors. Just as accelera- 
tors require larger and larger energies to probe smaller and smaller scales, 
more sensitive gravitational-wave detectors will be probing regions closer 
and closer to the horizon, potentially reaching Planck scales and beyond. 
What may be there, lurking? 
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Summary 


“The crushing of matter to infinite density by infinite tidal gravita- 
tion forces is a phenomenon with which one cannot live comfortably. 
From a purely philosophical standpoint it is difficult to believe that 
physical singularities are a fundamental and unavoidable feature of 
our universe [...] one is inclined to discard or modify that theory 
rather than accept the suggestion that the singularity actually occurs 
in nature.” 


Kip Thorne, Relativistic Stellar Structure and Dynamics (1966) 


“No testimony is sufficient to establish a miracle, unless the testi- 
mony be of such a kind, that its falsehood would be more miraculous 
than the fact which it endeavors to establish.” 


David Hume, An Enquiry concerning Human Understanding (1748) 


1 Introduction 


The discovery of the electron and the known neutrality of matter led in 1904 
to J. J. Thomson’s “plum-pudding” atomic model. Data from new scattering 
experiments was soon found to be in tension with this model, which was even- 
tually superseeded by Rutherford’s, featuring an atomic nucleus. The point-like 
character of elementary particles opened up new questions. How to explain the 
apparent stability of the atom? How to handle the singular behavior of the 
electric field close to the source? What is the structure of elementary particles? 
Some of these questions were elucidated with quantum mechanics and quantum 
field theory. Invariably, the path to the answer led to the understanding of 
hitherto unknown phenomena. 

The history of elementary particles is a timeline of the understanding of 
the electromagnetic (EM) interaction, and is pegged to its characteristic 1/r? 
behavior (which necessarily implies that other structure has to exist on small 
scales within any sound theory). 

Arguably, the elementary particle of the gravitational interaction are black 
holes (BHs). Within General Relativity (GR), BHs are indivisible and are in 
fact the simplest macroscopic objects that one can conceive. The uniqueness 
theorems — establishing that the two-parameter Kerr family of BHs describes 
any vacuum, stationary and asymptotically flat, regular BH solution — have 
turned BHs into somewhat of a miracle elementary particle gl. 

Can BHs, as envisioned in vacuum GR, hold the same surprises that the 
electron and the hydrogen atom did when they started to be experimentally 
probed? Are there any parallels that can be useful guides? The BH interior 
is causally disconnected from the exterior by an event horizon. Unlike the 
classical description of atoms, the GR description of the BH exterior is self- 
consistent and free of pathologies. The “inverse-square law problem” — the GR 


counterpart of which is the appearance of pathological curvature singularities — 
is swept to inside the horizon and therefore harmless for the external world. 
There are good indications that classical BHs are perturbatively stable against 
small fluctuations, and attempts to produce naked singularities, starting from 
BH spacetimes, have failed. BHs are not only curious mathematical solutions 
to Einstein’s equations, but also their formation process is well understood. In 
fact, there is nothing spectacular with the presence or formation of an event 
horizon. The equivalence principle dictates that an infalling observer crossing 
this region (which, by definition, is a global concept) feels nothing extraordinary: 
in the case of macroscopic BHs all of the local physics at the horizon is rather 
unremarkable. 

Why, then, should one question the existence of BHs? [] There are a number 
of important reasons to do so. 


e Even if the BH exterior is pathology-free, the interior is not. The Kerr 
family of BHs harbors, singularities and closed timelike curves in its in- 
terior, and more generically it features a Cauchy horizon signaling the 
breakdown of predictability of the theory. In fact, the geometry describ- 
ing the interior of an astrophysical spinning BH is currently unknown. A 
possible resolution of this problem may require accounting for quantum 
effects. It is conceivable that these quantum effects are of no consequence 
whatsoever on physics outside the horizon. Nevertheless, it is conceiv- 
able as well that the resolution of such inconsistency leads to new physics 
that resolves singularities and does away with horizons, at least in the 
way we understand them currently. Such possibility is not too dissimilar 
from what happened with the atomic model after the advent of quantum 
electrodynamics. 


e In a related vein, quantum effects around BHs are far from being under 
control, with the evaporation process itself leading to information loss. 
The resolution of such problems could include changing the endstate of 
collapse {4}/8| (perhaps as-yet-unknown physics can prevent the formation 
of horizons), or altering drastically the near-horizon region (9). The fact 
of the matter is that there is no tested nor fully satisfactory theory of 
quantum gravity, in much the same way that one did not have a quantum 
theory of point particles at the beginning of the 20th century. 


Tt is ironic that one asks this question more than a century after the first BH solution 
was derived. Even though Schwarzschild and Droste wrote down the first nontrivial regular, 
asymptotic flat, vacuum solution to the field equations already in 1916 Pi], several decades 
would elapse until such solutions became accepted. The dissension between Eddington and 
Chandrasekhar over gravitational collapse to BHs is famous — Eddington firmly believed that 
nature would find its way to prevent full collapse — and it took decades for the community 
to overcome individual prejudices. Progress in our understanding of gravitational collapse, 
along with breakthroughs in the understanding of singularities and horizons, contributed to 
change the status of BHs. Together with observations of phenomena so powerful that could 
only be explained via massive compact objects, the theoretical understanding of BHs turned 
them into undisputed kings of the cosmos. The irony lies in the fact that they quickly became 
the only acceptable solution. So much so, that currently an informal definition of a BH might 
well be “any dark, compact object with mass above three solar masses.” 


e It is tacitly assumed that quantum gravity effects become important only 
at the Planck scale. At such lengthscales \/ch/G, the Schwarzschild ra- 
dius is of the order of the Compton wavelength of the BH. In the orders 
of magnitude standing between the Planck scale and those accessible by 
current experiments many surprises can hide (to give but one example, 
extra dimensions would change the numerical value of the Planck scale). 


e Horizons are not only a rather generic prediction of GR, but their exis- 
tence is in fact necessary for the consistency of the theory at the classical 
level. This is the root of Penrose’s (weak) Cosmic Censorship Conjec- 
ture, which remains one of the most urgent open problems in fundamental 
physics. In this sense, the statement that there is a horizon in any space- 
time harboring a singularity in its interior is such a remarkable claim, that 
(in an informal description of Hume’s statement above) it requires similar 
remarkable evidence. This is the question we will entertain here, is there 
any observational evidence for the existence of BHs? 


e It is in the nature of science that paradigms have to be constantly ques- 
tioned and subjected to experimental and observational scrutiny. Most 
specially because if the answer turns out to be that BHs do not exist, the 
consequences are so extreme and profound, that it is worth all the possible 
burden of actually testing it. Within the coming years we will finally be in 
the position of performing unprecedented tests on the nature of compact 
dark objects, the potential pay-off of which is enormous. As we will argue, 
the question is not just whether the strong-field gravity region near com- 
pact objects is consistent with the Kerr geometry, but rather to quantify 
the limits of observations in testing the event horizon. 


e If we don’t entertain alternatives, we won’t find them. We need to know 
how to model alternatives to understand how consistent the data is with 
accepted paradigms. 


e Finally, there are sociological issues, sometimes invoked to claim that one 
should invest taxpayers money in sure-fire science. The reality is that 
taxpayers, all of us, want to know. More frequently than not, the thrill 
lies in the intellectual possibilities rather than on merely checking known 
facts. This is, after all, why pure science is relevant. Thus, it is a waste of 
resources not to use (built and already paid for) detectors to investigate 
issues that lie at the heart of fundamental questions. 


Known physics all but exclude BH alternatives. Nonetheless, the Standard 
Model of fundamental interactions is not sufficient to describe the cosmos — at 
least on the largest scales — and also leaves all of the fundamental questions 
regarding BHs open. It may be wise to keep an open mind. As we will argue, 
from the point of view of observers merely testing the spacetime structure with- 
out a priori theoretical bias, the evidence for horizons is nonexistent. In fact, 
the question to be asked is not if there is an event horizon in the spacetime, but 
how close to it do experiments or observations go. 


2 Setting the stage: escape cone, photospheres, 
quasinormal modes, and tidal effects 


? 


“Alas, I abhor informality.’ 
That Mitchell and Webb Look, Episode 2 


Our framework will be that of GR. Let us start by focusing on spherical sym- 
metry for simplicity. Birkhoff’s theorem guarantees that any vacuum spacetime 
is described by the Schwarzschild geometry, which in standard coordinates reads 


ds? = —fdt? + f-tdr? + r7d0?, (1) 


where M is the total mass of the spacetime (we use geometrical G = c = 1 
units, except if otherwise stated) and 


f=1- 5. (2) 


The existence of a null hypersurface at r = 2M is the defining feature of a 
BH. This hypersurface is called the event horizon of a Schwarzschild BH. In 
the following, we will compare the properties of the latter with those of an 
ultracompact object with a effective surface at 


ro = 2M(1+6), (3) 


having often in mind the case where e < 1 (for instance, for Planckian cor- 
rections of an astrophysical BH, e ~ 10740). Although the above definition is 
coordinate-dependent, the proper distance between the surface and r, scales like 

zry ar ft? ~ et/?. Most of the results discussed below show a dependence on 
log €, making the distinction irrelevant. There are objects for which the effective 
surface is ill-defined, because the matter fields are smooth everywhere. In such 


cases, ro can conservatively be taken to be the geometric center of the object. 


2.1 Geodesics 


The geodesic motion of timelike or null particles in the geometry can be de- 
scribed with the help of two conserved quantities, the specific energy E = f t and 
angular momentum L = r?¢%, where a dot stands for a derivative with respect 
to proper time pom]. The radial motion can be computed via a normalization 
condition, 

:3 2 L? = p2 

t =E IE rd) SB Veco 5 (4) 


r 


where 6; = 1,0 for timelike or null geodesics, respectively. The null limit can 
be approached letting E, L — oo and re-scaling all quantities appropriately. 
At large distances, r > M, the motion of matter resembles that in Newton’s 


theory. For instance, circular orbits, defined by 7 = 7 = 0, exist and are stable at 
large distance: small fluctuations in the motion restore the body to its original 
position. However, circular trajectories are stable only when r > 6M, and 
unstable for smaller radii. The r = 6M surface defines the innermost stable 
circular orbit (ISCO), and has an important role in controlling the inner part 
of the accretion flow onto compact objects. 


2.2 Escape trajectories 
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Figure 1: Left: Critical escape trajectories of radiation in the Schwarzschild 
geometry. A locally static observer (located at r = robs) emits photons isotropi- 
cally, but those emitted within the colored conical sectors will not reach infinity. 
Right: Coordinate roundtrip time of photons as a function of the emission angle 
W > Pesce and for e < 1. 


There are other GR effects on the motion of particles around compact ob- 
jects, beyond the appearance of ISCOs in spacetimes. One of them concerns the 
behavior of light rays and how they are bent by the spacetime curvature. There 
exists a critical value of the angular momentum L = KME for a light ray to 
be able to escape to infinity. By requiring that a light ray emitted at a given 
point will not find turning points in its motion, Eq. yields Kese = 3V3 (10). 
Suppose that the light ray is emitted by a locally static observer at r = rọ. In 
the local rest frame, the velocity components of the photon are 


MK 
go Z a/ fos quer — 1— K2M2 fo , (5) 


ro TO 


where fo = f(ro) = 1—2M/ro. With this, one can easily compute the escape 
angle, sin Yese = 3M V3 fo/ro. In other words, the solid angle for escape is 


27M2(rp — 2M) ro —2M 
SC = ies 2 —onar ? 
AQese = 20 (: y 1 3 a (6) 


where the last step is valid for e < 1. For angles larger than these, the light 
ray falls back and either hits the surface of the object, if there is one, or will 
be absorbed by the horizon. The escape angle is depicted in Fig. [1] for different 
emission points rọ. The rays that are not able to escape reach a maximum 
coordinate distance, 

A fo M? ) 


r2 sin? y 
This result is accurate away from Wesc, whereas for Y —> Yese the photon ap- 
proaches the photosphere (r = 3M) discussed in the next section. The coordi- 
nate time that it takes for photons that travel initially outward, but eventually 
turn back and hit the surface of the object, is shown in Fig. [1] as a function of 
the locally measured angle wv, and is of order ~ M for most of the angles wv, for 


e <1. We find a closed form expression away from Wet, which describes well 
the full range (see Fig. D, 


Tmax ~ 2M (1 + (7) 


troundtrip ~ 8M log (cot (w/2)) £ (8) 


When averaging over w, the coordinate roundtrip time is then 32M Cat/r ~ 
9.33M, for any «<1, where “Cat” is Catalan’s constant. 


2.3 Photospheres, ECOs and ClePhOs 
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Figure 2: Schematic classification of compact objects according to their com- 
pactness (in a logarithmic scale). Objects in the same category have similar 
dynamical properties on a timescale 7. 


Another truly relativistic feature is the existence of circular null geodesics, 
i.e., of circular motion for high-frequency electromagnetic or gravitational waves 
(GWs). In the Schwarzschild geometry, this can happen only at r = 3M. This 
location defines a surface called the photosphere, or, in the context of equatorial 
slices, a light ring. The photosphere has a number of interesting properties, and 


is useful to understand certain features of spacetimes. For example, it controls 
how BHs look like when illuminated by accretion disks or stars, thus defining 
their so-called “shadow”. Imaging these shadows for the supermassive dark 
source SgrA* is the main goal of the Event Horizon Telescope (13}[14]. The 
photosphere also has a bearing on the spacetime response to any type of high- 
frequency waves, and therefore describes how high-frequency gravitational waves 
linger close to the horizon. 

At the photosphere, V! = —2E?/(3M?) < 0. Thus, circular null geodesics 


geo 


are unstable: a displacement 6 of null particles grows exponentially 


At MEON oe eee 
S(t) ~ be", A= ETT (9) 


A geodesic description anticipates that light or GWs may persist at or close to 
the photosphere on timescales 33M ~ 5M. Because the geodesic calculation 
is local, these conclusions hold irrespectively of the spacetime being vacuum 
all the way to the horizon or not. An ultracompact object with surface at 
ro = 2M(1 + €), with € < 1, would feature exactly the same geodesics and 
properties close to its photosphere, provided that on timescales ~ 15M the null 
geodesics did not have time to bounce off the surface. We are requiring three 
e-fold times for the instability to dissipate more than 99.7% of the energy of the 
initial pulse. This amounts to requiring that 


E€ [Í cerit ~ 0.0165. (10) 


Thus, the horizon plays no special role in the response of high frequency waves, 
nor could it: it takes an infinite (coordinate) time for a light ray to reach the 
horizon. The above threshold on € is a natural sifter between two classes of 
compact, dark objects. For objects characterized by e€ 2 0.0165, light or GWs 
can make the roundtrip from the photosphere to the object’s surface and back, 
before dissipation of the photosphere modes occurs. For objects satisfying (LO), 
the waves trapped at the photosphere relax away by the time that the waves 
from the surface hit it back. 

We can thus use the properties of the ISCO and photosphere to distinguish 
between different classes of models: an object is defined as compact if it features 
an ISCO, or in other words if its surface satisfies rọ < 6M. Accretion disks 
around compact objects of the same mass should have similar characteristics. 
It is customary to define Exotic Compact Objects (ECOs) as those horizonless 
compact objects more massive than a neutron star. Among the compact objects, 
some feature photospheres. These could be called ultracompact objects (UCOs). 
Finally, those objects which satisfy condition have a “clean” photosphere, 
and will be designated by ClePhOs. The early-time dynamics of ClePhOs is 
expected to be the same as that of BHs. At late times, ClePhOs should display 
unique signatures of their surface. The zoo of compact objects is summarized 


in Fig. 


?This diagram can be misleading in special cases. Some objects with a surface at finite € can 


2.4 Quasinormal modes 


Low-frequency gravitational or EM waves are not well described by null par- 
ticles. In the linearized regime, massless waves are all described by a master 
partial differential equation of the form 


,OU(t,r) O?W(t,r) 
Or Ot? 


+ff 


V(r)U(t,r) = S(t,r). (11) 


The source term S(t, r) contains information about the cause of the disturbance 
W(t,r). The information about the angular dependence of the wave is encoded 
in the way the separation was achieved, and involves an expansion in scalar, 
vector or tensor harmonics for different spins s of the perturbation. These 
angular functions are labeled by an integer l > |s|, and the effective potential is 


Wi +1) 2, 2M 
vas (EPa) (12) 
with s = 0, +1, +2 for scalar, vector or axial tensor modes. The s = +2 equation 


does not describe completely gravitational perturbations, the potential for polar 
gravitational perturbations is more complicated [16}{18). 

It is hard to extract general information from the PDE (11), in particular 
because its solutions depend on the source term and initial conditions. We 
can gain some insight by studying the source-free equation in Fourier space. By 
defining the Fourier transform through W(t, r) = Tz J e~**th(w, r)dw, one gets 
the following ODE 

ay ?-V)y=0 13 


where r is implicitly written in terms of the new “tortoise” coordinate 
+2M1 ( $ 1) (14) 
z=r og | —— — ; 
8 OM 


such that z(r) diverges logarithmically near the horizon. In terms of z, Eq. 
is equivalent to the time-independent Schrödinger equation in one dimension 
and it reduces to the wave equation governing a string when M = l = 0. To 
understand a string of length L with fixed ends, one imposes Dirichlet boundary 
conditions and gets an eigenvalue problem for w. The boundary conditions can 
only be satisfied for a discrete set of normal frequencies, w = F (n = 1,2,...). 
The corresponding wavefunctions are called normal modes and form a basis onto 
which one can expand any configuration of the system. The frequency is purely 


real because the associated problem is conservative. 


nevertheless behave as ClePhOs. For example, constant density stars near the Buchdahl limit 
behave as ClePhOs for axial modes: these do not couple to the fluid and travel unimpeded to 
the center of the star. Their travel time can be considerably long, hence these stars act, for 
all purposes, as ClePhOs for these fluctuations. A specific example is discussed in footnote 
3. Nevertheless, it should be clear that these are contrived examples which are not meant to 
describe quantum corrections. 


10 


If one is dealing with a BH spacetime, the appropriate conditions (required 
by causality) correspond to having waves traveling outward to spatial infinity 
(W ~ e? as z — oc) and inwards to the horizon (Y ~ e~*** as z + —oo) [see 
Fig. B]. Due to backscattering off the effective potential (12), the eigenvalues 
w are not known in closed form, but they can be computed numerically 
(18). The fundamental | = 2 mode (the lowest dynamical multipole in GR) of 
gravitational perturbations reads 


Mw = M(wpr + iwr) © 0.373672 — 10.0889623 . (15) 


Remarkably, the entire spectrum is the same for both the axial or the polar 
gravitational sector (17). The frequencies are complex and are therefore called 
quasinormal frequencies. Their imaginary component describes the decay in 
time of fluctuations on a timescale T = 1/|wy|, and hints at the stability of the 
geometry. Unlike the case of a string with fixed end, we are now dealing with an 
open system: waves can travel to infinity or down the horizon and therefore it is 
physically sensible that any fluctuation damps down. The corresponding modes 
are the quasinormal modes (QNMs), which in general do not form a complete 
set. 

Boundary conditions play a crucial role in the structure of the QNM spec- 
trum. If a reflective surface is placed at ro = 2M(1+ €) Z 2M, where (say) 
Dirichlet boundary conditions have to be imposed, the spectrum changes con- 
siderably. For e = 10~° the fundamental mode of the eigenspectrum reads 


Mwpolar œ 0.13377 — 12.8385 x 1077 (16) 
Mwaxial œ 0.13109 — 12.3758 x 1077. (17) 


Not only are the two gravitational sectors no longer isospectral but, more im- 
portantly, the perturbations have smaller frequency and are much longer lived, 
since a decay channel (the horizon) has disappeared. 


More generically, for Dirichlet or Neumann boundary conditions at r = ro, 
the QNMs in the e > 0 limit read 


M 
Mwr ~ x~—(pr—5)~ | loge}, (18) 
2|z0| 
M = 
Mw, ~ -fig Mun) ~ —|log e| 7+3 , (19) 


where zo = z(ro) ~ 2M loge, p is an odd (even) integer for Dirichlet (Neumann) 
boundary conditions, ô is the phase of the wave reflected at r = ro, and Bis = 


2 
a ei The above scaling can be understood in terms of modes 
trapped between the peak of the potential at r ~ 3M and the “hard 
surface” at r = ro [see Fig. [B]. Low-frequency waves are almost trapped 
by the potential, so their frequency scales as the size of the cavity (in tortoise 
coordinates), wr ~ 1/zo, just like the normal modes of a string. The (small) 


imaginary part is given by waves which tunnel through the potential and reach 
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Figure 3: Typical effective potential for perturbations of a Schwarzschild BH 
(top panel) and of an horizonless compact object (bottom panel). In the BH 
case, QNMs are waves which are outgoing at infinity (z + +00) and ingoing at 
the horizon (z —> —oo), whereas the presence of a potential well (provided either 
by a reflective surface, a centrifugal barrier at the center, or by the geometry) 
supports quasi-trapped, long-lived modes. 


infinity. The tunneling probability can be computed analytically in the small- 
frequency regime and scales as |A|? ~ (Mwpr)”t? « 1 (23). After a time 
t, a wave trapped inside a box of size zg is reflected N = t/zo times, and 
its amplitude reduces to A(t) = Ao (1 — |A|?)* ~ Ao (1-t|Al?/z0). Since, 
A(t) ~ Age !#7!# ~ Ag(1 — |wz|t) in this limit, we immediately obtain 


wr ~ 1/20, wry LA|? /zo ~ wet? . (20) 


This scaling agrees with exact numerical results and is valid for any l and any 
type of perturbation. 

Clearly, a perfectly reflecting surface is an idealization. In certain models, 
only low-frequency waves are reflected, whereas higher-frequency waves probe 
the internal structure of the specific object 2723]. In general, the location of 
the effective surface and its properties (e.g., its reflectivity) can depend on the 
energy scale of the process under consideration. 
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2.5 Quasinormal modes, photospheres, and echoes 


The effective potential Vgeo for geodesic motion reduces to that of wave propa- 
gation V in the high-frequency, high-angular momentum (i.e., eikonal) regime. 
Thus, some properties of geodesic motion have a wave counterpart [1]. The 
instability of light rays along the null circular geodesic translates into some 
properties of waves around objects compact enough to feature a photosphere. 
A wave description needs to satisfy “quantization conditions.” Since GWs are 
quadrupolar in nature, the lowest mode of vibration should satisfy 


co 5 2 
Mw = 2% = —- ~ 0.3849. 21 

R i 3v3 ( ) 
In addition the mode is damped, as we showed, on timescales 3/3M. Overall 
then, the geodesic analysis predicts 


Mw ~ 0.3849 — 70.19245 . (22) 


This crude estimate, valid in principle only for high-frequency waves, matches 
rather well even the fundamental mode of a Schwarzschild BH, Eq. (15). 

Nevertheless, QNMs frequencies can be defined for any dissipative system, 
not only for compact objects or BHs. Thus, the association with photospheres 
has limits. Such an analogy is nonetheless enlightening in the context of objects 
so compact that they have photospheres and resemble Schwarzschild deep into 
the geometry, in a way that condition is satisfied [232429]. 

For a BH, it becomes clear that the excitation of the spacetime modes hap- 
pens at the photosphere [15]. The vibrations excited there travel outwards to 
possible observers or down the event horizon. The structure of GW signals at 
late times is therefore expected to be relatively simple. This is shown in Fig. |4 
which refers to the scattering of a Gaussian pulse of axial quadrupolar modes off 
a BH. The pulse crosses the photosphere, and excites its modes. The ringdown 
signal, a fraction of which travels to outside observers, is to a very good level 
described by the lowest QNMs. The other fraction of the signal generated at 
the photosphere travels downwards and into the horizon. It dies off and has no 
effect on observables at large distances. 

Contrast the previous description with the dynamical response of ultracom- 
pact objects for which condition is satisfied (i.e., a ClePhO) [cf. Fig. [4]. 
The initial description of the photosphere modes still holds, by causality. Thus, 
up to timescales of the order |zọ| ~ —M loge (the roundtrip time of radiation 
between the photosphere and ro) the signal is identical to that of BHs (23) [24]. 
At later times, however, the pulse traveling inwards is bound to see the object 
and be reflected either at its surface or at its center. In fact, this pulse is semi- 
trapped between the object and the light ring. Upon each interaction with the 
light ring, a fraction exits to outside observers, giving rise to a series of echoes 
of ever-decreasing amplitude. From Eqs. (18)—(19), repeated reflections occur 
in a characteristic echo delay time 


Techo ~ 4M| log el . (23) 
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However, the main burst is typically generated at the photosphere and has 
therefore a frequency content of the same order as the BH QNMs (15). The 
initial signal is of high frequency and a substantial component is able to cross 
the potential barrier. Thus, asymptotic observers see a series of echoes whose 
amplitude is getting smaller and whose frequency content is also going down Ë] 
At very late times, the QNM analysis is valid, and the signal is described by 
Eqs. (18)—(19). 

A formal description of the the frequency content of the sequence of echoes, 
and how they can be obtained with the help of the BH response, can be found 
in Ref. [26]. A method to reconstruct the potential of the echo source, after the 
modes have been extracted was proposed in Ref. [33]. 


2.6 The role of the spin 


When spherical symmetry is broken by angular momentum J, Birkhoff’s the- 
orem does not apply, so the exterior metric of a spinning object is not unique 
and technically much more difficult to computeļff] Furthermore, BH uniqueness 
theorems do not apply in the presence of matter. In particular, a stationary 
object is not necessarily axisymmetric, unlike the Kerr solution which uniquely 
describes a BH in stationary configurations. Even restricting to axisymmet- 
ric configurations, the linear response of a spinning object depends also on the 
azimuthal number m (such that |m| < 1), which characterizes the angular de- 
pendence of the perturbations and yields a Zeeman-like splitting of the QNM 
frequencies. The relation between null geodesics and BH QNM s in the eikonal 
limit is more involved but conceptually similar to the static case [86]. 

Assuming that the exterior spacetime approaches the Kerr geometry when 
e€ — 0, the linear response of a ClePhO can be studied with a toy model 
in which the Kerr horizon is replaced by a reflective surface at ro = r+(1+ €), 
where rz = M(1 + \/1— x?) are the locations of the horizons and y = J/M? 
is the dimensionless spin. In this case, Eq. is modified by the replacement 
WR — wR — MQ on the left-hand side (where Q is the angular velocity of 
the object when € —> 0), whereas Eq. is replaced by 


) beier (wR—mO), (24) 


s [ 2M? 
Mwy œ~ Bi (aoe 


~ [zo Ty 77 


3 An interesting, nontrivial example was worked out years ago in a different context [Bo]. 
The example concerns axial GWs emitted during the scatter of point particles around ultra- 
compact, constant-density stars. The GW signal shows a series of visible, well-space echoes 
after the main burst of radiation, which are associated to the quasi-trapped s-modes of ul- 
tracompact stars BI. It turns out that all the results of that work fit very well our descrip- 
tion: axial modes do not couple to the fluid and travel free to the geometrical center of the 
star, which is therefore the effective surface in this particular case. The time delay of the 
echoes in Fig. 1 of Ref. is very well described by the GW’s roundtrip time to the center, 
Techo ~ 222e71/2M, where ro = 2M(1 + €) is the radius of the star and rọ = 2M is the 


Buchdahl’s limit zi 
4Nonetheless, there are indications that all multipole moments of the external spacetime 


approach those of a Kerr BH as « > 0 [34][35]. In this limit, it is natural to expect that the 
exterior geometry is extremely close to Kerr, unless some discontinuity occurs in the BH limit. 
It would be extremely important to confirm or confute this conjecture. 
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Figure 4: Ringdown waveform for a BH (dashed black curve) compared to a 
ClePhO (solid red curve) with a reflective surface at ro = 2M(1 + e) with 
e = 1071}. We considered | = 2 axial gravitational perturbations and a Gaussian 
wavepacket w(r,0) = 0, w(r, 0) = e7 (7m)? /0° (with zm = 9M and o = 6M) 
as initial condition. Note that each subsequent echo has a smaller frequency 
content and that the damping of subsequent echoes is much larger than the 
late-time QNM prediction (e~“%!* with wr M ~ 4 x 1071? for these parameters). 
Data available online {19]. 


where now zo ~ M[1+(1—yx?)~!/?] loge. The result shows how angular mo- 
mentum can bring about substantial qualitative changes. The spacetime is un- 
stable for wg(wr— mQ) < 0 (i.e., in the superradiant regime z ), on a timescale 
Tinst = 1/wz. This phenomenon is called ergoregion instability . In the 
e — 0 limit and for sufficiently large spin, wr ~ MQ and wy; ~ | log e|~*. Notice 
that such description is only valid at small frequencies, and therefore becomes 
increasingly less accurate at large spins and away from the instability threshold 
(where wgr,wr ~ 0). For large spin the result is more complex and can be found 
in Ref. [43]. Furthermore, in the superradiant regime the “damping” factor, 
wr/wr > 0, so that, at very late times (when the pulse frequency content is 
indeed described by these formulas), the amplitude of the QNMs increases due 
to the instability. Such increase is anyway small, for example = x 4x 1076 
when e = 0.001, / = 2 and y = 0.7. 
In the spinning case, the echo delay time reads 


Tecno ~ 2M[1+ (1 — x?)7/?] loge, (25) 


which corresponds to the period of the corotating mode, Techo ~ (wr —mQ)7!. 
Note that, as we explained earlier, the signal can only be considered as a se- 
ries of well-defined pulses at early stages. In this stage, the pulse still contains a 
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substantial amount of high-frequency components which are not superradiantly 
amplified. Thus, amplification will occur at late times; the early-time evolution 
of the pulse generated at the photosphere is more complex. 


2.7 BHs in binaries: tidal heating and tidal deformability 


The properties of an event horizon have also important consequences for the 
dynamics of binary systems containing a BH. As previously discussed, a spinning 
BH absorbs radiation of frequency w > mQ, but amplifies radiation of smaller 
frequency. In this respect, BHs are dissipative systems which behave just like 
a Newtonian viscous fluid [aljas]. Dissipation gives rise to various interesting 
effects in a binary system — such as tidal heating [46147], tidal acceleration, and 
tidal locking — the Earth-Moon system being no exception due to the friction of 
the oceans with the crust. 

For low-frequency circular binaries, the energy flux associated to tidal heat- 
ing at the horizon corresponds to the rate of change of the BH mass [48|[49], 


5 

M = By x CE Or -2), (26) 

where Qg « 1/M is the orbital angular velocity and the (positive) prefactor 

depends on the masses and spins of the two bodies. Thus, tidal heating is 

stronger for highly spinning bodies relative to the nonspinning by a factor ~ 
0/QK > 1. 

The energy flux leads to a potentially observable phase shift of GWs 
emitted during the inspiral. Thus, it might be argued that an ECO binar 
can be distinguished from a BH binary, because Ey = 0 for the former 50). 
However, the trapping of radiation in ClePhOs can efficiently mimic the effect 
of a horizon [50]. In order for absorption to affect the orbital motion, it is 
necessary that the time radiation takes to reach the companion, Taa, be much 
longer than the radiation-reaction time scale due to heating, Ter ~ E/Eu, 
where E ~ —}M(MQx)?/ is the binding energy of the binary (assuming equal 
masses). For BHs, Tjaa + co because of time dilation, so that the condition 
Traa >> Tre is always satisfied. For ClePhOs, Taq is of the order of the GW 
echo delay time, Eq. (23), and therefore increases logarithmically as € — 0. 
Thus, an effective tidal heating might occur even in the absence of an event 
horizon if the object is sufficient compact. The critical value of e increases 
strongly as a function of the spin. For orbital radii larger than the ISCO, the 
condition Traq >> Trr requires € < 10788 for y < 0.8, and therefore even Planck 
corrections at the horizon scale are not sufficient to mimic tidal heating. This 
is not necessarily true for highly spinning objects, for example Taa > TRR at 
the ISCO requires e < 10716 for y ~ 0.9. 

Finally, the nature of the inspiralling objects is also encoded in the way they 
respond when acted upon by the external gravitational field of their companion — 
through their tidal Love numbers [51]. An intriguing result in classical GR is 
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the fact that the tidal Love numbers of a BH are precisely zero 54|. On 
the other hand, those of ECOs are small but finite [55]. In particular, the 
tidal Love numbers of ClePhOs vanish logarithmically in the BH limit. Thus, 
any measurement of the tidal Love number k translates into an estimate of the 
distance of the ECO surface from its Schwarzschild radius, 


enve F, (27) 


Owing to the above exponential dependence, a measurement of the tidal Love 
number of an ECO at the level of k ~ O(107?) already probes Planck distances 
away from the gravitational radius r, [55]. This level of accuracy is in principle 
within reach future GW detectors 5o 
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3 Beyond vacuum black holes 


“Mumbo Jumbo is a noun and is the name of a grotesque idol said to 
have been worshipped by some tribes. In its figurative sense, Mumbo 
Jumbo is an object of senseless veneration or a meaningless ritual.” 

Concise Oxford English Dictionary 


To entertain the possibility that dark, massive, compact objects are not BHs, 
requires one to discuss some outstanding issues. One can take two different 
stands on this topic: 


i. a pragmatic approach of testing the spacetime close to compact, dark 
objects, irrespective of their nature, by devising model-independent ob- 
servations that yield unambiguous answers. 


ii. a less ambitious and more theoretically-driven approach, which starts by 
constructing objects that are very compact, yet horizonless, within some 
framework. It proceeds to study their formation mechanisms and stability 
properties; then discard solutions which either do not form or are unsta- 
ble on relatively short timescales; finally, understand the observational 
imprints of the remaining objects, and how they differ from BHs’. 


In practice, when dealing with outstanding problems where our ignorance 
is extreme, both approaches should be used simultaneously. Indeed, using con- 
crete models can sometimes be a useful guide to learn about broad, model- 
independent signatures. As it will hopefully become clear, one could design 
exotic horizonless models which mimic all observational properties of a BH with 
arbitrary accuracy. While the statement “BHs exist in our Universe” is funda- 
mentally unfalsifiabld?| alternatives can be ruled out or confirmed to exist with 
a single observation, just like Popper’s black swans. 


3.1 Are there alternatives? 


A nonexhaustive summary of possible objects which could mimic BHs is shown 
in Table Stars made of constant density fluids are perhaps the first known 
example of compact configurations, the literature on the subject is too long to 
list. In GR, isotropic static spheres made of ordinary fluid satisfy the Buchdahl 
limit on their compactness, 2M/ro < 8/9 [82], and can thus never be a ClePhO. 
Interestingly, polytrope UCO stars always have superluminal sound speed [57]. 

Compact solutions can be built with fundamental, massive bosonic fields. 
The quest for self-gravitating bosonic configurations started in an attempt to 


5 Arguably, this is true for any statement about the existence of a physical entity in our 
Universe, but this does not prevent us to accumulate more and more evidence to support a 
given model or theory, nor to rule out competitors. 
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understand if gravity could produce “solitons” through the nonlinearities of 
the field equations. Such configurations are broadly referred to as boson stars. 
Boson stars have attracted considerable interest as light scalars are predicted to 
occur in different scenarios, and ultralight scalars can explain the dark matter 
puzzle (136). When fastly spinning, such stars can be extremely compact. There 
seem to be no studies on the classification of such configurations (there are 
solutions known to display photospheres, but it is unknown whether they can 
be as compact as ClePhOs). 

The other objects in the list require either unknown matter or large quantum 
effects. For example, a specific model of a wormhole was discussed carefully in 
the context of “BH foils” [102]. In the nonspinning case, the geometry of this 
model is described by 


ds? = —(f +A) dt? 4 “ + dQ? (28) 


The constant À is assumed to be extremely small, A ~ e-M*/ lp where fp is 
the Planck length. In the context of GR, this solution requires matter violating 
the null energy condition. Other objects, such as fuzzballs, were introduced 
as a microstate description of BHs in string theory. In this setup the indi- 
vidual microstates are horizonless and the horizon arises as a coarse-grained 
description of the microstate geometries. Phenomena such as Hawking radia- 
tion can be recovered, in some instances, from classical instabilities (111}/137]. 
“Gravitational-vacuum stars” or gravastars are ultracompact configura- 
tions supported by a negative pressure, which might arise as an hydrodynamical 
description of one-loop QFT effects in curved spacetime, so they do not neces- 
sarily require exotic new physics [138]. In these models, the Buchdahl limit is 
evaded because the internal effective fluid is anisotropic and also violates the 
energy conditions (139). Gravastars have been recently generalized to include 
anti-de Sitter cores, in what was termed AdS bubbles, and which may allow for 
holographic descriptions [o8]. 

For most objects listed in Table[I}which are inspired by quantum-gravity cor- 
rections, the changes in the geometry occur deep down in the strong field regime. 
Some of these models — including also black stars [124], superspinars or 
collapsed polymers — predict that large quantum backreaction should 
affect the horizon geometry even for macroscopic objects. In these models, our 
parameter € is naturally of the order ~ O(€p/M) € (10739, 10746) for masses 
in the range M € (10,108)Mo. The 2 — 2-hole model predicts even more com- 
pact objects, with € ~ (€p/M)? € (10778,1079?) |118|. In all these cases, both 
quantum-gravity or microscopic corrections at the horizon scale select ClePhOs 
as well-motivated alternatives to BHs. 

Despite a number of supporting arguments — some of which urgent and well 
founded — it is important to highlight that there is no horizonless ClePhO for 
which we know sufficiently well the physics at the moment. 
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3.2 Formation and evolution 


Although supported by sound arguments, the vast majority of the alternatives 
to BHs are, at best, incompletely described. Precise calculations (and often 
even a rigorous framework) incorporating the necessary physics are missing. 

One notable exception to our ignorance are boson stars. These configurations 
are known to arise, generically, out of the gravitational collapse of massive 
scalars. Their interaction and mergers can be studied by evolving the Einstein- 
Klein-Gordon system, and there is evidence that accretion of less massive boson 
stars makes them grow and cluster around the configuration of maximum mass. 
In fact, boson stars have efficient gravitational cooling mechanisms that allow 
them to avoid collapse to BHs and remain very compact after interactions 
a. 

The remaining objects listed in Table |1| were built in a phenomenological 
way or they arise as solutions of Einstein equations coupled to exotic matter 
fields. For example, models of quantum-corrected objects do not include all 
the (supposedly large) local or nonlocal quantum effects that might prevent 
collapse from occurring. In the absence of a complete knowledge of the missing 
physics, it is unlikely that a ClePhO forms out of the merger of two ClePhOs. 
These objects are so compact that at merger they will be probably engulfed by 
a common apparent horizon. The end product is, most likely, a BH. On the 
other hand, if large quantum effects do occur, they would probably act on short 
timescales to prevent apparent horizon formation. In some models, Planck- 
scale dynamics naturally leads to abrupt changes close to the would-be horizon, 
without fine tuning (118). Likewise, in the presence of (exotic) matter or if GR 
is classically modified at the horizon scale, Birkhoff’s theorem no longer holds, 
and a star-like object might be a more natural outcome than a BH. In summary, 
with the exception of boson stars, we do not know how or if UCOs and ClePhOs 
form in realistic collapse or merger scenarios. 
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Model Taxonomy Formation Stability EM signatures GWs 
Fluid stars UCOs V v 
Anisotropic stars ClePhOs x V v V 
Boson stars & UCOs, (ClePhOs?) v vV vV v 
Gravastars COs — ClePhOs x V v ~ 
AdS bubbles UCOs — ClePhOs x vV ~ x 
Wormholes ClePhOs x V v ~ 
99}{103 
Fuzzballs ClePhOs x x ~ 
Blfe] (but see [110}f113}) (but see [23|[24]/114)) 
Superspinars COs — ClePhOs x v x ~ 
[115] 116 (but see (117) 
2 — 2 holes ClePhOs x x x ~ 
118 (but see [118}) (but see [118}) 
Collapsed ClePhOs x x ~ 
polymers (but see [119}[121]) 
Quantum bounces /| ECO — ClePhOs x x ~ 
black stars (but see [123][126]) 125 
Quantum stars* UCOs — ClePhOs x x x x 
127/128 
Fire-walls* ePhOs x x x ~ 
129H131 132 


Table 1: Catalogue of some proposed horizonless compact objects. A v tick 
means that the topic was addressed. With the exception of boson stars, how- 
ever, most of the properties are not fully understood yet. An asterisk * stands 
for the fact that these objects are BHs, but could have phenomenology similar 
to the other compact objects in the list. This table does not include models of 
quantum-corrected BHs (e.g., [133}{135)), even when the latter predict large cor- 
rections near the horizon, because also in this case the object is a BH (although 
different from Kerr). 
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3.3 On the stability problem 


“There is nothing stable in the world; uproar’s your only music.” 
John Keats, Letter to George and Thomas Keats, Jan 13 (1818) 


Appealing solutions are only realistic if they form and remain as long-term 
stable solutions of the theory. In other words, solutions have to be stable when 
slightly perturbed or they would not be observed. There are strong indications 
that the exterior Kerr spacetime is stable, although a rigorous proof is still miss- 
ing [t40]. There are good reasons to believe that some — if not all — horizonless 
compact solutions are linearly or nonlinearly unstable. 

Some studies of linearized fluctuations of ultracompact objects are given in 
Table We will not discuss specific models, but we would like to highlight 
two general results. Linearized gravitational fluctuations of any nonspinning 
UCO are extremely long-lived and decay no faster than logarithmically 
mjs]. Indeed, such perturbations can be again understood in terms of 
modes quasi-trapped within the potential barrier shown in Fig.|3} they require 
a photosphere but are absent in the BH case. For a ClePhO, these modes are 
very well approximated by Eqs. (18)—(19). The long damping time of these 
modes has led to the conjecture that any UCO is nonlinearly unstable and may 
evolve through a Dyson-Chandrasekhar-Fermi type of mechanism B64]. The 
endstate is unknown, and most likely depends on the equation of the state of the 
particular UCO: some objects may fragment and evolve past the UCO region 
into less compact configurations, via mass ejection, whereas other UCOs may 
be forced into gravitational collapse to BHs. 

The above mechanism is supposed to be active for any spherically symmetric 
UCO, and also on spinning solutions. However, it is nonlinear in nature. On the 
other hand, UCOs (and especially ClePhOs) can develop negative-energy regions 
once spinning. In such a case, there is a well-known linear instability, dubbed 
as ergoregion instability. The latter affects any horizonless geometry with an 
ergoregion and is deeply connected to superradiance [22]. 

In summary, there is good evidence that UCOs are linearly or nonlinearly 
unstable. Unfortunately, the effect of viscosity is practically unknown [56], and 
so are the timescales involved in putative dissipation mechanisms that might 
quench this instability. On the other hand, even unstable solutions are relevant 
if the timescale is larger than any other dynamical scale in the problem [} The 
nonlinear mechanism at work for spherically symmetric spacetimes presumably 
acts on very long timescales only; a model problem predicts an exponential 
dependence on the size of the initial perturbation (143]. At least in some portion 
of the parameter space the ergoregion-instability timescale is very large 
[79]. From Eq. (24), we can estimate the timescale of the instability of a spinning 


6 After all, “we are all unstable”. 
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ClePhO, 


+1- x) (= — =) lwa(r4 -r Ch) (29) 


2Bis Fa WR- mQ 


As previously discussed, a spinning ClePhO is (superradiantly) unstable only 
when Q > wr/m. For example, for l = m = s = 2 and y = 0.7, the above 
formula yields 


1 1 
T= —~—|loge| 
WI 


M 
r € (5,1) (aar) yr when e€ (10~*°,10-°”). (30) 


Generically, the ergoregion instability acts on timescales which are parametri- 
cally longer than the dynamical timescale, ~ M, of the object. Given such long 
timescales, it is likely that the instability can be efficiently quenched by some 
dissipation mechanism of nongravitational nature, although this effect would be 
model-dependent [89]. Furthermore, it is possible that the endstate of the in- 
stability is simply an ECO spinning at the superradiant threshold, Q = wg/m. 
In such case, the instability would only rule out highly-spinning ECOs in a 
certain compactness range. Likewise, EM or GW observations indicating sta- 
tistical prevalence of slowly-spinning compact objects, across the entire mass 
range, could be an indication for the absence of an horizon. 

Finally, there are indications that instabilities of UCOs are merely the equiv- 
alent of Hawking radiation for these geometries, and that therefore there might 
be a smooth transition in the emission properties when approaching the BH 


limit (102}[111]. 
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4 Have exotic compact objects been already ruled 
out by electromagnetic observations? 


“And here - ah, now, this really is something a little recherché.” 
Sherlock Holmes, The Musgrave Ritual, Sir Arthur Conan Doyle 


Before the GW revolution, questions about the observational evidence for 
BHs were initially posed in the context of EM observations. A concise overview 
of the arguments is presented in Refs. [144145], concluding that “it is funda- 
mentally impossible to give an observational proof for the existence of a BH 
horizon.” Surprisingly however, there have been recent claims that 
observations of some objects (most notably the compact radio source Sgr A*) all 
but rule out other candidates, including those whose surface is a Planck-distance 
away from the horizon (i.e., € ~ 10745 for Sgr A*, thus qualifying as ClePhOs 
according to our classification). The argument can, roughly, be described as 
follows: 


i. There are systems with a very low accretion rate and very low luminosity, 
for which the central object is very dark. 


ii. The coordinate time that a particle in the accretion disk takes to reach 
the surface scales as Tiraver ~ M| log e|. Even for Planck-sized € this time 
is very small. 


iii. Assumption: because the above timescale is so small, a thermodynamic 
and dynamic equilibrium must be established between the accretion disk 
and the central object, on relatively short timescales. 


iv. Assumption: Because there is equilibrium, the central object must be 
returning in EM radiation most of the energy that it is taking in from the 
disk. 


v. Because we see no such radiation, the premise that there is a surface must 
be wrong. 


The arguments above fail to take into account important physics: 


Timescale to reach equilibrium. Strong lensing can prevent equilibrium 
from being achieved around ultracompact objects in sufficiently short time. 
Consider matter in the accretion disk falling in, releasing scattered radiation 
on the surface of the object which is then observed by our detectors. Sup- 
pose, for the sake of the argument, that once hitting the surface, it is scattered 
isotropically. Then, as discussed in Section [2.2] only a fraction ~ e is able to 
escape during the first interaction with the star, cf. Eq. (6). The majority of 
the radiation, will fall back onto the surface after a time troundtrip ~ 9.3M given 
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by the average of Eq. u Suppose one injects, instantaneously, an energy 6M 
onto the object. Then, after a time T, the energy emitted to infinity during 
N =T/trounatrip interactions reads 


AE ~ [1-(1-6)*] ôM. (31) 


Clearly, AE => 6M as N —> oo since all energy will eventually escape to infinity. 
However, AE ~ eNOM if eN <1, i.e. when 


z M tHubble 
107168 32 
E (saa) ( T ) ? (a 


where we have been extremely conservative in normalizing T by the age of the 
Universe, tHubble © 101° yr. The luminosity of ClePhOs satisfying the above 


property is roughly 
. 17 € ôM 
Bed a) (St) : oy 


When e < 10716, this quantity is minute and it is impossible to achieve equilib- 
rium on any meaningful timescale. Indeed, any injected energy 6M is eventually 
radiated in a timescale longer than tHubble Whenever Eq. is satisfied. Thus, 
gravastars, fuzzballs, and other objects inspired by quantum effects at the hori- 
zon scale (e < 10~*°) are not ruled out as exotic supermassive objects. 
Cascade into other channels. Even if the object were returning all of 
the incoming radiation on a sufficiently short timescale, a sizeable fraction of 
this energy could be in channels other than EM. Let us estimate the center- 
of-mass energy in collisions between matter being accreted and the compact 
object. The infalling matter is assumed to be on radial free fall with four- 
velocity vA) = (E/f,—y/ E? — f,0,0). For particles at the surface of the object 
via) = (/f,0,0,0). When the particles collide, their CM energy reads 149), 


= [Se Mov 2B 
Ecm = mov2 1— Juv (ya) oe a ; (34) 


with mo the mass of the colliding particles (assumed to be equal). Therefore, 
the Lorentz factor of the colliding particles, in the CM frame, is y >> 3 x 10° for 
e€ X 10777, and y ~ 10!! for Planck-scale modifications of supermassive objects. 
Therefore, Ecm is orders of magnitude above the maximum energy achievable 
at the LHC. At these CM energies, it is possible that new degrees of freedom 
are excited (some of which might even be the unknown physics that forms these 
compact objects in the first place). It is well possible that if new fields exist 
they carry a substantial fraction of the luminosity (while going undetected in our 


TOne might wonder if the trapped radiation bouncing back and forth the surface of the 
object might not interact with the accretion disk. As we showed in Section [2.2] this does not 
happen, as the motion of trapped photons is confined to within the photosphere. 
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telescopes $] Even without advocating new physics beyond the 10 TeV scale, 
extrapolation of known hadronic interactions to large energies suggests that 
about 20% of the collision energy goes into neutrinos, whose total energy is a 
sizeable fraction of that of the photons emitted in the process (151). 

In summary, we believe that some observations of supermassive objects most 
likely are compatible with ClePhOs (in particular they do not rule out gravastars 
nor quantum corrections at the horizon scale); these same observations are 
probably not compatible with the remaining of the UCO parameter space. 

No argument whatsoever can be made to rule out exotic compact objects of 
arbitrary compactness, simply because the transition to the BH limit should be 
continuous. Therefore, any argument suggesting that exotic compact objects 
are excluded for any € has to be taken with great care. 


8Incidentally, gravitational radiation makes only negligible contribution to the total flux, 
even for such large CM energies. The ratio of EM to gravitational radiation generated during 
these processes was estimated to be [150] 


Egem logy q 
Egrav Y- mo í 


(35) 


with q/mo the charge-to-mass ratio of the colliding particles. Thus, even at such large CM 
energies, EM radiation is vastly larger than gravitational radiation, because q/mọo ~ 101° for 
protons. 
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5 Testing the nature of BHs with GWs 


“Tt is well known that the Kerr solution provides the unique solu- 
tion for stationary BHs in the universe. But a confirmation of the 
metric of the Kerr spacetime (or some aspect of it) cannot even be 
contemplated in the foreseeable future.” 


S. Chandrasekhar, The Karl Schwarzschild Lecture, 
Astronomische Gesellschaft, Hamburg (September 18, 1986) 


Testing the nature of dark, compact objects with EM observations is plagued 
with several difficulties. Some of these are tied to the incoherent nature of the 
EM radiation in astrophysics, and the amount of modeling and uncertainties 
associated to such emission. Other problems are connected to the absorption 
by the interstellar medium. As discussed in the previous section, testing quan- 
tum or microscopic corrections at the horizon scale with EM probes is very 
challenging. 

The historical detection of GWs by aLIGO opens up the exciting pos- 
sibility of testing gravity in extreme regimes with unprecedented accuracy 
[132)[153}156]. GWs are generated by coherent motion of massive sources, and 
are therefore subjected to less modeling uncertainties (they depend on far fewer 
parameters) relative to EM probes. The most luminous GWs come from very 
dense sources, but they also interact very feebly with matter, thus providing 
the cleanest picture of the cosmos, complementary to that given by telescopes 
and particle detectors. 

Compact binaries are the preferred sources for GW detectors (and in fact 
were the source of the first detected events [1524157] ). The GW signal from com- 
pact binaries is naturally divided in three stages, corresponding to the different 
cycles in the evolution driven by GW emission [158}{160]: the inspiral stage, 
corresponding to large separations and well approximated by post-Newtonian 
theory; the merger phase when the two objects coalesce and which can only be 
described accurately through numerical simulations; and finally, the ringdown 
phase when the merger end-product relaxes to a stationary, equilibrium solution 
of the field equations [16}{160)/161]. 

All three stages provide independent, unique tests of gravity and of compact 
GW sources. The GW signal in the two events reported so far is consistent 
with them being generated by BH binaries, and with the endstate being a BH. 
To which level, and how, are alternatives consistent with current and future 
observations? This question was recently addressed by several authors and may 
be divided into two different schemes (see Table [2). 


5.1 Smoking guns: echoes and resonances 


For binaries composed of ClePhOs, the GWs generated during inspiral and 
merger is expected to be very similar to those by a corresponding BH binary with 
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BH ECO ClePhO 
Ss 
3 GW echoes x v (only UCOs) V (Techo ~ M| log el) 
oh Modified prompt ringdown x v x 
= Extra modes x vV v 
la 
m Multipolar structure (2PN) 6M, = 6S; = 0 6M, #0, 6S; #0 6M, ~ 0, 6S; ~ 0 
A=] Tidal heating (2.5 — 4PN) v x x 
5 Tidal Love number (5PN) k=0 k < O(kns) k ~ [loge]~! 
S Resonances x 80\|96| wM ~ [log e]—+ 


Table 2: Summary of the properties of exotic compact objects that affect the 
GW signal relative to BHs in the ringdown phase and in the inspiral phase at 
a given post-Newtonian (PN) order. In the multipolar structure, 6M; and ôS; 
are the corrections to the mass and current_multipole moments relative to the 
no-hair relation for a Kerr BH in isolation (162), Mı + iS; = M'*! (ix)!. Tidal 
absorption at the horizon enters at 2.5PN (4PN) order for spinning (nonspin- 
ning) BHs. The typical tidal Love numbers of a neutron star are denoted by 
kys. 


the same mass and spin. The arguments were detailed in Section [2] However, 
clear distinctive features appear due to the absence of a horizon. These features 
are: 


e The appearance of late-time echoes in the waveforms (23][24|/29][40]. After 
the merger, ClePhOs give rise to two different ringdown signals: the first 
stage is dominated by the photosphere modes, and is indistinguishable 
from a BH ringdown. However, after a time Techo ~ M|loge| (see Sec- 
tion 2-4), the waves trapped in the “photosphere+ClePhO” cavity start 
leaking out as echoes of the main burst. This is a smoking-gun for new 
physics, potentially reaching microscopic or even Planckian corrections at 


the horizon scale (23}[24). 

Strategies to dig GW signals containing “echoes” out of noise are not 
fully under control, first efforts are underway [163}/164). The ability to 
detect such signals depends on how much energy is converted from the 
main burst into echoes (i.e., on the relative amplitude between the first 
echo and the prompt ringdown signal in Fig. |4). Define the ratio of en- 
ergies to be Yecho- Then, the signal-to-noise ratio p necessary for echoes 
to be detectable separately from the main burst, for a detection thresh- 


old of p = 8, is -i > ——82__.. In such case the first echo 
P Pprompt ringdown ~< Vecho(%) 


is detectable with a simple ringdown template (an exponentially damped 
sinusoid [165}). Space-based detectors will see prompt ringdown events 
with very large p |166]. For Yeho = 20%, we estimate that the planned 
space mission LISA will see at least one event per year, even for 
the most pessimistic population synthesis models used to estimate the 
rates (166). The proposed Einstein telescope or Voyager-like [169] 
third-generation Earth-based detectors will also be able to distinguish 
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ClePhOs from BHs with such simple-minded searches. The event rates 
for LIGO are smaller, and more sophisticated searches need to be im- 
plemented. Preliminary analysis of GW data using the entire echoing 
sequence was reported recently (38){170|/171]. More detailed characteriza- 
tion of the echo waveform (e.g., using Green’s function techniques and 
accounting for spin effects |40|) is necessary to reduce the systematics in 
data analysis. In several models of ClePhOs, yYecho is large enough to pro- 
duce effects that — also owing to the logarithmic dependence of Techo for 
all models known so far — are within reach of near-future GW detectors, 
even if the ClePhO corrections occurs at the Planck scale [23][24]. This is 
a truly remarkably prospect. As the sensitivity of GW detectors increases, 
the absence of echoes might be used to rule out many quantum-gravity 
models completely, and to push tests of gravity closer and closer to the 
horizon scale, as now routinely done for other cornerstones of GR, e.g. in 


tests of the equivalence principle [1554/172]. 


e Model-dependent fluid modes are also excited. Due to redshift effects, 
these will presumably play a subdominant role in the GW signal. 


e Certain models of ECOs arise naturally in effective theories with extra 
gravitational degrees of freedom . In such case, the detection of extra 
polarizations (as achievable in the future with a global network of inter- 
ferometers) might provide evidence for new physics at the horizon scale. 


e Resonant mode excitation during inspiral. The echoes at late times are 
just vibrations of the ClePhO. These vibrations have relatively low fre- 
quency and can, in principle, also be excited during the inspiral process 
itself, leading to resonances in the motion as a further clear-cut signal of 


new physics ; 


e Spin-mass distribution of compact objects skewed towards low spin, across 
the mass range. The development of the ergoregion instability depletes 
angular momentum from spinning ClePhOs, independently of their mass. 
Although the effectiveness of such process is not fully understood, it would 
lead to slowly-spinning objects as a final state. 


e It is also possible that, at variance with the BH case, the merger of two 
ClePhOs can be followed by a burst of EM radiation associated with the 
presence of high-density matter during the collision. Although such emis- 
sion might be strongly redshifted, searches for EM counterparts of candi- 
date BH mergers can provide another distinctive signature of new physics 
at the horizon scale. 


5.2 Precision physics: QNMs, tidal deformability, heat- 
ing, and multipolar structure 


The GW astronomy era will also gradually open the door to precision physics, 
for which smoking signs may not be necessary to test new physics. 
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e If the product of the merger is an ECO but not a ClePhO, it will simply 
vibrate differently from a BH. Thus, precise measurements of the ring- 
down frequencies and damping times allow one to test whether or not 
the object is a BH {174]/175]. Such tests are in principle feasible for wide 


classes of objects, including boson stars {80 , gravastars [95[971139], 
121125]. 


wormholes , or other quantum-corrected objects 


e The merger phase can also provide information about the nature of the 
coalescing objects. ECOs which are not sufficiently compact will likely 
display a merger phase resembling that of a neutron-star merger rather 
than a BH merger (24][178}. The situation for ClePhOs is unclear since no 
simulations of a full coalescence are available. 


e The absence of a horizon affects the way in which the inspiral stage pro- 
ceeds. In particular, three different features may play a role, all of which 
can be used to test the BH-nature of the objects. 


— No tidal heating for ECOs. Horizons absorb incoming high frequency 
radiation, and serve as sinks or amplifiers for low-frequency radiation 
able to tunnel in (see Section [2.7p. UCOs and ClePhOs, on the 
other hand, are not expected to absorb any significant amount of 
GWs. Thus, a “null-hypothesis” test consists on using the phase of 
GWs to measure absorption or amplification at the surface of the 
objects (50). LISA-type GW detectors will place stringent tests 
on this property, potentially reaching Planck scales near the horizon 
and beyond [50]. 


— Nonzero tidal Love numbers for ECOs. In a binary, the gravitational 
pull of one object deforms its companion, inducing a quadrupole 
moment proportional to the tidal field. The tidal deformability is 
encoded in the Love numbers, and the consequent modification of 
the dynamics can be directly translated into GW phase evolution at 
higher-order in the post-Newtonian expansion . It turns out that 
the tidal Love numbers of a BH are zero [52}[53][180] 183|, allowing 
again for null tests. By devising these tests, existing and upcoming 
detectors can rule out or strongly constraint boson stars 
184| or even generic ClePhOs [50455]. 


— Different multipolar structure. Spinning objects in a binary are also 
expected to possess multipole moments that differ from those of the 
Kerr geometry. This property is true at least for UCOs and for less 
compact objects. The impact of the multipolar structure on the GW 
phase will allow one to estimate and constrain possible deviations 
from the Kerr geometry {185}. 
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6 Summary 


“The important thing is not to stop questioning. Curiosity has its 
own reason for existing.” 


Albert Einstein, From the memoirs of William Miller, an editor, 
quoted in Life magazine, May 2, 1955; Expanded, p. 281 


Thomson’s atomic model was carefully constructed, and tested theoretically 
for inconsistencieq)| Rutherford’s incursion into scattering of œ particles were 
not meant to disprove the model, they were aimed at testing its accuracy. Ac- 
cording to Marsden, “...it was one of those ’hunches’ that perhaps some effect 
might be observed, and that in any case that neighbouring territory of this 
Tom Tiddler’s ground might be explored by reconnaissance. Rutherford was 
ever ready to meet the unexpected and exploit it, where favourable, but he also 
knew when to stop on such excursions” [187]. 

All the current observational evidence gathered around massive, compact 
and dark objects is compatible with the BH hypothesis. There is no equally- 
well-motivated alternative, satisfying known laws of physics and composed of 
ordinary matter, which is compatible with those same observations. Despite 
this, there are long-standing problems associated with horizons and singulari- 
ties, which hint at some inconsistency between classical gravity and quantum 
mechanics at the scale of the horizon. Furthermore, since the majority of the 
matter in our Universe is unknown, we cannot exclude that dark compact ob- 
jects made of exotic matter are lurking in the cosmos. Even, and if, these issues 
are resolved within classical physics, some outstanding questions remain. What 
evidence do we have that compact objects are BHs? How deep into the poten- 
tial can observations probe, and up to where are they still compatible with the 
current paradigm? At one hand’s reach, we have the possibility to dig deeper 
and deeper into the strong-field region of compact objects, for free. As the sen- 
sitivity of EM and GW detectors increases, so will our ability to probe regions 
of ever increasing redshift. Perhaps the strong-field region of gravity holds the 
same surprises that the strong-field EM region did? 


It is amusing that some alternative models — which we now know to be closer to Ruther- 
ford’s model — were studied, shown to be unstable and therefore ruled out [186]. 
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